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Damping of magnetization dynamics in a ferromagnetic metal, arising from spin-orbit coupling, 
is usually characterised by the Gilbert parameter a. Recent calculations of this quantity, using 
a formula due to Kambersky, find that it is infinite for a perfect crystal owing to an intraband 
scattering term which is of third order in the spin-orbit parameter This surprising result conflicts 
with recent work by Costa and Muniz who study damping numerically by direct calculation of 
the dynamical transverse susceptibility in the presence of spin-orbit coupling. We resolve this 
inconsistency by following the approach of Costa and Muniz for a slightly simplified model where 
it is possible to calculate a analytically. We show that to second order in ^ one retrieves the 
Kambersky result for a, but to higher order one does not obtain any divergent intraband terms. 
The present work goes beyond that of Costa and Muniz by pointing out the necessity of including 
the effect of long-range Coulomb interaction in calculating damping for large A direct derivation 
of the Kambersky formula is given which shows clearly the restriction of its validity to second order 
in f so that no intraband scattering terms appear. This restriction has an important effect on the 
damping over a substantial range of impurity content and temperature. The experimental situation 
is discussed. 


I. INTRODUCTION 

Magnetization dynamics in a ferromagnetic metal is central to the field of spintronics with its many applications. 
Damping is an essential feature of magnetization dynamics and is usually treated phenomenologically by means 
of a Gilbert term in the Landau-Lifshitz-Gilbert equation Qi- For a system with spin-rotational invariance the 
uniform precession mode of the magnetization in a uniform external magnetic field is undamped and the fundamental 
origin of damping in ferromagnetic resonance is spin-orbit coupling (SOG). Early investigations of the effect include 
those of Kambersky M and Korenman and Prange Kambersky’s torque-correlation formula for the Gilbert 
damping parameter a has been used by several groups [Ml, some of whom have given alternative derivations. 
However the restricted validity of this formula, as discussed below, has not been stressed. In this torque-correlation 
model contributions to a of both intraband and interband electronic transitions are usually considered. The theory 
is basically developed for a pure metal with the effect of defects and/or phonons introduced as phenomenological 
broadening of the one-electron states. Assuming that the electron scattering-rate increases with temperature T due 
to electron-phonon scattering the intraband and interband transitions are found to play a dominant role in low and 
high T regimes, respectively. The intraband (interband) term is predicted to decrease(increase) with increasing T and 
to be proportional to where ^ is the SOG parameter. Accordingly a is expected to achieve a minimum at an 

intermediate T. This is seen experimentally in Ni and hep Co [l^ but not in Fe [T^ and FePt [l^. The ^ dependence 
of a is well-established at high T but there seems to be no experimental observation of the predicted ^ 

behaviour at low T. The interband term in Kambersky’s theory can be given a very simple interpretation in terms of 
second-order perturbation theory Q . A quite different phenomenological approach, not applicable in some unspecified 
low scattering-rate regime, has been adopted to try and find a physical interpretation of the intraband term [^, . 

No acceptable theoretical treatment of damping in this low scattering regime is available because the intraband term 
of Kambersky’s theory diverges to infinity in the zero-scattering limit of a pure metal with translational symmetry at 
T=0 |9l.Illl Il3l|. We consider it essential to understand the pure metal limit before introducing impurity and phonon 
scattering in a proper way. 

Costa and Muniz [l9j | recently studied damping numerically in this limit by direct calculation of the dynamical 
spin susceptibility in the presence of SOG within the random phase approximation (RPA). They determine a from 
the linewidth of the uniform (wave-vector q = 0) spin-wave mode which appears as a resonance in the transverse 
susceptibility. One of the main objects of this paper is to establish some degree of consistency between the work of 
Kambersky and that of Costa and Muniz. We follow the approach of the latter authors for a slightly simplified model 
where it is possible to calculate a analytically. We show that to second order in ^ one retrieves the Kambersky result, 
but to higher order no intraband terms occur, which removes the problem of divergent a. To confirm this point, in 
Appendix we derive the Kambersky formula directly in a way that makes clear its restriction to second order in 
^ to which order the divergent terms in a arising from intraband transitions do not appear. This throws open the 
interpretation of the minimum observed in the temperature dependence of a for Ni and Co. 

At this point we may mention an alternative theoretical approach to the calculation of Gilbert damping using 
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scattering theory [2^ [^. Starikov et al [2l| find that, for Nii_xFex alloys at T=0, a becomes large near the pure 
metal limits x=0,l. They attribute this to the Kambersky intraband contribution although no formal correspondence 
is made between the two approaches. 

The work of Costa and Muniz [l^ follows an earlier paper where it is shown that SOC has the effect of 
coupling the transverse spin susceptibility to the longitudinal spin susceptibility and the charge response. It is known 
that a proper calculation of these last two quantities in a ferromagnet must take account of long-range Coulomb 
interactions (23 - [T7| . The essential role of these interactions is to ensure conservation of particle number. Costa et 
al [l^, [13 do not consider such interactions but we show here that this neglect is not serious for calculating a with 
sufficiently small SOC. However in the wider framework of this paper, where mixed charge-spin response is also readily 
studied, long-range interactions are expected to sometimes play a role. They also come into play, even to second order 
in when inversion symmetry is broken. 

In section ni we establish the structure of spin-density response theory in the presence of SOC by means of exact 
spin-density functional theory in the static limit . In section IIIII we introduce a spatial Fourier transform and an 
approximation to the dynamical response is obtained by introducing the frequency dependence of the non-interacting 
susceptibilities. The theory then has the same structure as in the RPA. Section H^is devoted to obtaining an explicit 
expression for the transverse susceptibility in terms of the non-interacting susceptibilities. Expressions for these, in 
the presence of SOC, are obtained within the tight-binding approximation in section [V] In section IVll we consider 
the damping of the resonance in the q=0 transverse susceptibility and show how the present approach leads to the 
Kambersky formula for the Gilbert damping parameter a where this is valid, namely to second order in the SOC 
parameter We do not give an explicit formula for a beyond this order but it is clear that no intraband terms appear. 
In section I VIII some experimental aspects are discussed with suggestions for future work. The main conclusions are 
summarized in section rVIIIl 


II. SPIN-DENSITY FUNCTIONAL THEORY WITH SPIN-ORBIT COUPLING 

The Kohn-Sham equation takes the form 

(v) + = en(t>ua{r) ( 1 ) 

with the spin index a =t, \. corresponding to quantization along the direction of the ground-state magnetization in a 
ferromagnet. This may be written in 2 x 2 matrix form with eigenvectors ■ The density matrix is defined 

in terms of the spin components (/'rio-(r) of the one-electron orbitals by 

— 'y '^ 4^na{j^')4^na' ^n) (2) 

n 

where 0[x) is the unit step function and /rg is the chemical potential. The electron density is given by 

(7 ncT 

and the effective potential in ([1]) is 

v;^//(r) = Waa' (r) -k 6aa' J (fr'p{r')v{r - r') -k (r) (4) 

where Waa'i'r) is the external potential due to the crystal lattice and any magnetic fields and u(r) = e^/|r| is the 
Coulomb potential. The exchange-correlation potential v^^,{r) is defined as 6Exc/Sncra'{^), a functional derivative of 
the exchange-correlation energy E^c- The term in ([IJ is the SOC energy. A small external perturbation dw^^>, 
for example due to a magnetic field, changes the effective potential to -k , giving rise to new orbitals and 
hence to a change in density matrix . The equation 

6naa>{r) = (5) 

where H is the volume of the sample, defines a non-interacting response function and the full response function x 
is defined by 

Sucrcr'ir) = J d^r'xa,7’TT'(r,r')SwrT'(r'}. 


( 6 ) 
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From (|4]) 


= Sw„^a[irl) +J2 j d^r2[viri - (r 2 ) 


( 7 ) 


and we may write 


^^(72(72 ('' 2 ) 


52 ^ 


^^(7272 (^ 2 )^^Cri(T^ (^ 1 ) 


T^CTicr'0-2(72 *' 2 )- 


( 8 ) 


Combining ® - ® we find the following integral equation for the spin-density response function Xaa'TT'{i^, r'): 

Xo(7'rr'(r,r') = x°^,^^,(r,r') - (ff)"! ^ ^ 

(T 1 (T (T2^2 

X727^TT' (I ’2 I I" ) ■ 

(9) 


dVi / dV 2 X°„,„^^,(r,ri)[u(ri-r 2 )^, 


717^^7272 


K. 


(Tl(T^(T2<72 


(l’l:r 2 )] 


This equation is a slight generalisation of that given by Williams and von Barth [ 2 ^ . In the static limit it is formally 
exact although the exchange-correlation energy E^c is of course not known exactly. In the next section we generalise 
the equation to the dynamical case approximately by introducing the frequency dependence of the non-interacting 
response functions also take a spatial Fourier transform. In the case where SOC is absent the result is directly 

compared with results obtained using the RPA. 


III. DYNAMICAL SUSCEPTIBILITIES IN THE PRESENCE OF SPIN-ORBIT COUPLING AND 

LONG-RANGE GOULOMB INTERACTION 

In general the response functions x(r, r') are not functions of r — r' and a Fourier representation of ® for a 
spatially periodic system involves an infinite number of reciprocal lattice vectors. There are two cases where this 
complication is avoided. The first is a homogeneous electron gas and the second is in a tight-binding approximation 
with a restricted atomic basis. We may then introduce Fourier transforms of the form x(r) = X]q or 

x(q) = / d^rxir)e~'^'’^ and write ® as 

Xo-ct'tt'(Qj w) = X(T(t'tt'(Q! ~ ^ ^ ^ ^ Xo-(T'(Ti(tJ (*ll ^)^lo-io-2cr2 (*l)Xo-20-27r'(Q: (Id) 

(T2<^2 


where we have also introduced the uj dependence of x as indicated at the end of the last section. Here U(q) is an 
ordinary Fourier transform, without a factor so that 

hcTi (7^(72(72 (q) r^(q)^(7i(7^^(72(72 E (H) 

where w(q) = Aire'^/q^ is the usual Fourier transform of the Coulomb interaction and the second term is independent 
of q since AT is a short-range spatial interaction. In the gas case it is proportional to a delta-function 5(r — r') in the 
local-density approximation (LDA) and in tight-binding it can be taken as an on-site interaction. In both cases 
K may be expressed in terms of a parameter U as 

d^7i7'^7272 17[Jcri(7^ ^(72(72 ^(71 (72 T (7^ ^CT2 (7^ ^(7^ (72 ] (^2) 

where a = 4 ,, t for a ='|',in the tight-binding case this form of K corresponds to a simple form of interaction 
which leads to a rigid exchange splitting of the bands ( [ 2 ^, [l^). This is only appropriate for transition metals in a 
model with d bands only, hybridization with s and p bands being neglected. We adopt this model in order to obtain 
transparent analytic results as far as possible. Although not as realistic as ’’first-principles” models of the electronic 
structure it has been used, even with some quantitative success, in treating the related problem of magnetocrystalline 
anisotropy in Co/Pd structures as well as pure metals [s^. In (fTOl) the response functions x are per unit volume 
in the gas case but, more conveniently, may be taken as per atom in the tight-binding case with v{q) modified to 
?;(q) = /(q^ila) where ila is the volume per atom. 
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To show how equations (nni) - (HU) are related to RPA we examine two examples in the absence of SOC. First 
consider the transverse susceptibility X 4 ,'t“|- 4 ,(q, w) which is more usually denoted by y_|_(q, w). Equation (fTUl) becomes 

xtm = x°tt4. “ (13) 

and, from CD and CD, Hut = -^UU = —U. Hence 

Xitti = X°n4.(l - ^ x°tu)”^ (14) 

which is just the RPA result of Izuyama et al [slj for a single-orbital Hubbard model and of Lowde and Windsor [slj 
for a five-orbital d-band model. Clearly in the absence of SOC the Coulomb interaction u(q) plays no part in the 
transverse susceptibility, as is well-known. A more interesting case is the longitudinal susceptibility denoted by Xmm 
in the work of Kim et al ( M, [ 13 ) and in . This involves only the response functions Xo-o-tt which we abbreviate 
to Xctt- In fact [23 


Xmm = xn + XU - XU - XU- (15) 

Without SOC xH takes the form XctHt and CD becomes 

Xcrr = Xo-Hr ~ Xcr Hrcr 2 Xo' 2 'r (15) 


with 14 £,2 = '^^(q) “ U 6 aa 2 ■ On solving the 2x2 matrix equation (fTOl) for x<tt, and using (ITSl) . we find the longitudinal 
susceptibility in the form 

_ _ X° + X° - - 2u(q)]x?X° _ 

Xmm j (-^0 X°)b(q) -U] + U[U - 2-(;(q)]x°X° 

which agrees with the RPA result that Kim et al ( [l^], [13) found for a single-orbital model. The Coulomb interaction 
v{q) is clearly important, particularly for the uniform susceptibility with q=0, where v U oo. It plays an essential 
role in enforcing particle conservation and hence in obtaining the correct result of Stoner theory. In view of the 
correspondence between our approach and the RPA method it seems likely that when SOC is included our procedure 
using equations (CD - CD should be almost equivalent to that of Costa and Muniz [l3 in the case of a model with 
d-bands only. However our inclusion of the long-range Coulomb interaction will modify the results. 


IV. AN EXPLICIT EXPRESSION FOR THE TRANSVERSE SUSCEPTIBILITY 


In this section we obtain an explicit expression for the transverse susceptibility Xtttt terms of the non-interacting 
response functions x°- We consider equation CD as an equation between 4x4 matrices where crtr' =j,t, tJ-, tt, 4-1 
labels the rows in that order and tt' labels columns similarly. The formal solution of (fTUl) is then 

X=(1 + X°H)-H°- (18) 


This expression could be used directly as the basis of a numerical investigation similar to that of Costa and Muniz. 
However we wish to show that the present approach leads to a Gilbert damping parameter a in agreement with the 
Kambersky formula, to second order in the SOC parameter ^ where Kambersky’s result is valid. This requires some 
quite considerable analytic development of (ITSl) . 

First we partition each matrix in (fTUl) into four 2x2 matrices. Thus from (fTTl) and (fTUl) 


with 


V = 






Ha 


fv — U V \ 

\ V v-u) ' 


(19) 


( 20 ) 


X = 


Also 


Xll Xl2 
X21 X22 


( 21 ) 
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and similarly for If we write 


dUD becomes 


where 


1 + = 


V X 21^1 


X? 2^2 \ 

1 + X 22 ^ 2 / 


A B 
C D 


/ 5-1 


/vO 

i All 

Xl2^ 


D-^ +D-^CS-^BD-'^) 

too 

too 

to 


S = A + BD-^C. 


( 22 ) 


(23) 


(24) 


The transverse susceptibility XittI which we are interested is the top right-hand element of xii so this is the 
quantity we wish to calculate. From 

Xii = ^-'(x?i-i?i^-'x2i) (25) 


and, from (l24ll and (l22|) . 


5=1+ X?il5 - X?2(l"2“' + X22)-'X2ll^l- 


(26) 


The elements of the 2x2 matrix S are calculated by straight-forward algebra and 

S'!! = 1 - t/x°tt4 + + x°u4)xtmXitn - (^ + Xtm)xunXmt 

+ X4,4.'t-|-)X't-|-'t-4,X4,f4,4, + + X't-|-t-|-)X4,4.'t-|-X4,'|-4.4.] 

where 


X = [v-U]/[U{U-2v)], Y = -v/[U{U -2v)] (28) 

and 

^ = (-’^ + Xtttt)("^ + X°4.4.4.) - (1^ + Xtt4.4.)(^ + Xutt) (29) 

The other three elements of S are given in Appendix |B] The transverse susceptibility is obtained from (EH) as 

Xim = [522 (x?i - - 5i2(x?i - BD-\°M/iSiiS 22 - 5i252i) (30) 

and 

B1^-'=X?2(1"2-'+X22)''- (31) 

A comparison of the fairly complex equation above for the transverse susceptibility with the simple well-known result 
(HI shows the extent of the new physics introduced by SOC. This is due to the coupling of the transverse susceptibility 
to the longitudinal susceptibility and the charge response, both of which involve the long-range Coulomb interaction. 

To proceed further it is necessary to specify the non-interacting response functions Xaa'aia' which occur throughout 
the equations above. 


V. THE NON-INTERACTING RESPONSE FUNCTIONS 

In the tight-binding approximation the one-electron basis functions are the Bloch functions 

|k^cr) = iV-l/2 gik-R, f-^2) 

3 

where j and ^ are the site and orbital indices, respectively, and N is the number of atoms in the crystal. The 
Hamiltonian in the Kohn-Sham equation now takes the form 

k/J,I/'CT 




( 33 ) 
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where T^i, corresponds to electron hopping and 

= -(a/2)(A + 6ex) (34) 

where u = 1,—1 for spin t,)- respectively. Here A = 2U{S^)/N where is the total spin angular momentum, in 
units of h, and the Zeeman splitting bex = where B^x: is the external magnetic field and is the Bohr 

magneton. The spin-orbit term ^ hj ■ Sj takes the second-quantized form 

= (^/2) ^ (35) 


where are matrix elements of the atomic orbital angular momentum operators ± iL^ in units 

of h. Within the basis of states ((32)l eigenstates of take the form 


and satisfy the equation 


fl<7 


H^ff\\in) =E^^\kn). 


Thus 


“"k/icr 




kn 


(36) 


(37) 


(38) 


where creates the eigenstate |kn). 

The non-interacting response function x^cr'aia' (*4’ ‘^) conveniently expressed as the Fourier transform of a retarded 
Green function by the Kubo formula 




"k+q/j, cr'" 


E ^ki-'fTi ^k+qi-'crj ))a 


(39) 


where the right-hand side is to be evaluated using the one-electron Hamiltonian Consequently, using ([38)l . we 

have 


xLvkt; (q> w) = E E + q)*<i.(k)a^Uk)*a^i^(k -f- q)((c k+qm'^kn; Ck„Ck-|-qm) ),, 




^ ^ ^ £^k+qm - -huEVf] 

kflu ran ^ 


(40) 


The last step uses the well-known form of the response function per atom for a non-interacting Fermi system (e.g. [1^) 
and 77 is a small positive constant which ultimately tends to zero. The occupation number /kn = F'(i?kn ~ Mo) where 
F is the Fermi function with chemical potential /io- Clearly for q = 0 the concept of intraband transitions (m = n), 
frequently introduced in discussions of the Kambersky formula, never arises for finite w since the difference of the 
Fermi functions in the numerator of (1401) is zero. Equation (1401) may be written in the form 


xLvi,Ti(q.w) = N 


-1 


kmn 


-Bmn(k, q)S^VHk, q)’' 


/k n fk+c 


Ek+qm F\z_n Eo -(- ITj 


(41) 


where 


■B™n(k,q) = ^a^^(k-hq)*a^^(k). 

M 


( 42 ) 
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VI. FERROMAGNETIC RESONANCE LINEWIDTH; THE KAMBERSKY FORMULA 


We now consider the damping of the ferromagnetic resonance in the <7 = 0 transverse susceptibility. The present 
approach, like the closely-related one of Costa and Muniz [T^, is valid for arbitrary strength of the SOC and can 
be used as the basis of numerical calculations, as performed by the latter authors. However it is important to show 
analytically that the present method leads to the Kambersky Q formula for the Gilbert damping parameter where 
this is valid, namely to second order in the SOC parameter This is the subject of this section. 

It is useful to consider first the case without SOC = 0). The eigenstates n of then have a definite spin and 
may be labelled na. It follows from (HOI) that X^aa'ai<j{ ^ ^crcrj^cr'CTi• Hence Xi 2 = 0 and, from (IHTl) . BD~^ = 0. Also, 
from Appendix B, S 12 = S '21 = 0. Thus, dSOD reduces to (ITTl) as it should. Considering x°.|-.|' 4 ,( 0 , w), given by (HOI) 
and m, we note that state m is pure j, spin, labelled by m and n is pure labelled by n f- Hence for ^ = 0 


from closure. Thus 


= ^(km|k/r)(k/r|kn) = Sjun 


kn 


/kn^ /k n\. 

Ewni - Ekn^ - hio + ir] 


and it follows from (IM)) that Eif^na- may be written as 

Euncr = Akn ~ (<^/2)(A -|- bex)- 

Hence we find from m that for ^ = 0 

= {2{S^)/N){bex -huj + iTj)-\ 


(43) 


(44) 


(45) 

(46) 


Thus, as ry —>• 0, i3x4.tt4.(*^!‘^) ^ sharp delta-function resonance at huj = b^x as expected. 

When SOC is included fvjj acquires an imaginary part that corresponds to damping. We now proceed to calculate 
this imaginary part to O(C^). To do this we can take ^ = 0 in the numerator of (l30t so that 

X4.tt4.(0)‘*^) = “ 5 'i25'2i/*S'22) (47) 

In fact S '12 and S '21 are both 0(^^) while S 22 is 0(1). Thus to obtain huj to 0(^^) we need only solve Sn = 0. 
Furthermore all response functions such as Xiftj,) with all but one spins in the same direction, are zero for ^ = 0 and 
need only be calculated to O(^) in (l?7l) . We show below that to this order they vanish, so that to 0(^^) the last term 
in Sii is zero and we only have to solve the equation 1 — = 0 for hio. This means that to second order in ^ 

the shift in resonance frequency and the damping do not depend on the long-range Coulomb interaction. 

To determine x°-|-'|'4,(0, w) to first order in ^ from (|40)) we notice that states n must be pure t spin, that is |kn) = 
|kri t)) while states m must be calculated using perturbation theory. The latter states may be written 

|kml) = |kTO t) ^E^^V ptr) (48) 

^ -C/kmt “ -^kptr 


_ £/kmi — -C^kptT 
pa ^ ^ 


where we have put and to first order in 


^tm “ icEE(' 


t * i 


/knt /kmt 


t * i 


/knt /kmt 


N ^ E - - feu + iry ^ E _ _ Scu + 

k/ii/ mn ( ( / I « 


(49) 


-) (50) 


with = (k/icr|kms), s = 1, 2, and = (kp,|kn) is independent of spin. Since atii/ ~ C we take = Om/i in 
the first term of dSOl)- Also aE^anfj. = by closure so that the first term of (15(11) vanishes since the difference 
of Fermi functions is zero. Only the second term of X^tl remains and this becomes, by use of (H51) . 


Xtni 


(kpt 4-)* * _ *_ /knt - /kml 

Akmt-Akpt + 

Kfiv mnp t' I I 


(51) 









Again using closure only terms with p = m = n survive and the matrix element of becomes 

(kn t I X! ^ {kn\L- |kn) = 0 (52) 

3 

due to the quenching of total orbital angular momentum L = Lj [s^. Thus, to first order in w), and 

similar response functions with one reversed spin, are zero. Hence we have only to solve 1 — Uxi-fn = 0 to obtain 
S(fia;) to 0(^^). Here we assume the system has spatial inversion symmetry without which the quenching of orbital 
angular momentum, as expressed by no longer pertains [s^ ■ We briefly discuss the consequences of a breakdown 
of inversion symmetry at the end of this section. 

On introducing the perturbed states (1481) and (|4^ we write dm in the form 




km n 


/k nl /kml 

-^kml -S^knl 

/knl /km2 


/kn2 /kml 


I oit 12 _ 

' E^,^ 2 -Eum-hio + irj 


' - i?krr 2 - So; + i77 

I rjtt 12 _ fkn2 ~ /km2 _, 

■' E^jr, 2 -Ekn 2 -hlU + i7^’' 


(53) 


Clearly and are of order B^^-^ is 0{^^) and is 0(1). We therefore neglect the term 

and, using (im and (im . we find 


-^rnlnl ^m 2 n 2 


^ (km|L |kn) 
2 T'knt A^km- 


kmf 


(54) 


The evaluation of |H^ 2 niP requires more care. It appears at first sight that to obtain this to 0(^^) we need to include 
second order terms in the perturbed eigenstates given by (H51) and (HHl) . However it turns out that these terms do not 
in fact contribute to \B~ 


_ 2 _j|^ to 0{^^) so we shall not consider them further. Then we find 

_ (km|L^|kn) (kw|L^|kp)(kp| L^|krt) 

™ 2"1 4 4 " {Ekrn - E^p) - E^p) 


(55) 


and hence to 0(^^) 


4 |(krn|L^|kp)p 2 l(kw|T^|kn)p 


I r„2nll ^ ^ (£;^_ _ ^^_)2 i ^ 


{Eun — ^^km)^ 


(56) 


The contribution of this quantity to y°.j..|.j^(0, w) in (l5^ may be written to 0{^^) as 


J_ IrH 

jy- / ^ I'^m 2 nl 


/knl /km2 


-( 1 - 


bex - fmj 


km 7 


-£'km2 -£^knl ^ex -£'km2 -S'knl ^ex H“ 


)■ 


(57) 


This is obtained by introducing the identity —Huj = —bex+{bex—Eu) in the relevant denominator in (I53L and expanding 
to first order in bex — but which turns out to be 0(4)- The remaining factors of this second term in (1571) may then 
be evaluated with ^ = 0, as at the beginning of this section, so that this term becomes {huj — bex)/{2U^{S^)). By 
combining equations (ESI), dSl and dm, and ignoring some real terms, we find that the equation 1 — U xj..|.|j^(0, w) = 0 
leads to the relation 


'^{tuj) = 7r^^/(2(S'^)) ^ [(/knt “ hrnl)\{km\L''\'kn)f5{El^rnl - Ekn^ - bex) 
kn m 

“l“(l/4)(/kn't' /kn4' /kmt /km^.) | l^^)l '^(-^km -S'kn ^ea:)] 

The Gilbert damping parameter a is given by ^{huj)/bex (e.g. [1^) and in (1551) we note that 

(/kn'I' /km4.)^(-^km4. T'kn'l' bex) — [7^(7?knt Mo) 7^(A^knt T ^ex A^o)]^(A'km4. T'kn'f' bex) 

— bexb{E\^ji^ /ro)^5(7?km'f' Mo) 


(58) 


(59) 


to Hrst order in bex at temperature T = 0. Similarly 


(/knCT fkma) b {E]^yji E\^xi bex) — ^ea,'^(A'kncr Mo)^(A'km(T Mo)- 


(60) 
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Thus from (l58)) 


a = ^eV(2(^")) ^ \{km\L^\kn)\HiEun^ - fio)6{E^mi - /io) 

kn™ 

+7^^^/(8(>S"')) ^ |(km|L \kn)\'^S{Ei,ric, - no)S{Ekrna - ^J'o) 

kn rn(T 

correct to 0(5^). We note that there is no contribution from intraband terms since (kn|L|kn) = 0. It is straight¬ 
forward to show that to O(C^) this is equivalent to the expression 

a = V(2(S'))EE l^nicr^ncr' {k)\^8{E^rna " ” Mo) (62) 

kn m gg ' 


where 


24mcr,n<7'(k) = C(kTOcr|[5' ,h®°]|knCT') 

and S~ is the total spin operator Sj with Sj = SJ — iS'j'. This may be written more concisely as 

a = 7t/{ 2{S^)) Y, \Amn(k)fSiEkm - Mo)<5(i;k„ - Mo) 

knm 


(63) 


(64) 


with 


24m„(k) = f(km|[S' ,/i®°]|kn) (65) 

and the understanding that the one-electron states km, kn are calculated in the absence of SOC. Equation (|M)) is the 
standard form of the Kambersky formula ( i, i) but in the literature SOC is invariably included in the calculation 
of the one-electron states. This means that the intraband terms with m = n no longer vanish. They involve the square 
of a delta-function and this problem is always addressed by invoking the effect of impurity and/or phonon scattering 
to replace the delta-functions by Lorentzians of width proportional to an inverse relaxation time parameter . Then 
as one approaches a perfect crystal (t —>■ oo) the intraband contribution to a tends to infinity. This behaviour is 
illustrated in many papers ( a, 0, El, m ). In fig. 1 of [ll it is shown clearly that a remains finite if one does 
not include SOC in calculating the one-electron states. However the effect of not including SOC is not confined to 
total removal of the intraband contribution. The remaining interband contribution is increased considerably in the 
low scattering rate regime, by almost an order of magnitude in the case of Fe. This makes a almost independent of 
scattering rate in Fe which may relate to its observed temperature independence [ll . The corresponding effect in 
Co is insufficient to produce the increase of a at low scattering rate inferred from its temperature dependence. The 
non-inclusion of SOC in calculating the one-electron states used in the Kambersky formula clearly makes a major 
qualitative and quantitative change in the results. This occurs as soon as intraband terms become dominant in 
calculations where they are included. For Fe, Co and Ni this corresponds to impurity content and temperature such 
that the scattering rate 1/r due to defects and/or phonons is less that about 10^"^sec“^ ( BUI)- Typically these 
metals at room temperature find themselves well into the high scattering-rate regime where the damping rate can 
be reliably estimated from the Kambersky interband term, with or without SOC included in the band structure [s^ . 
The physics at room temperature is not particularly interesting. One needs to lower the temperature into the low 
scattering-rate regime where intraband termsOf t hey exist, will dominate and lead to an anomalous ^ dependence of 
the damping on spin-orbit parameter ^ ( 0, [I3l)- The origin of this behaviour is explained in 0, [T^. It arises 

in the k sum of (|64|) from a striplike region on the Fermi surface around a line where two different energy bands cross 
each other in the absence of SOC. The strip width is proportional to ^ , or more precisely |^|. Since H„„(k) is of 
order ^ the contribution of intraband terms in (1641) is proportional to Thus the intraband terms lead to terms 
in a which diverge in the limit 0 and are non-analytic functions of The calculation of a in this section can 

be extended to higher powers of ^ than the second. No intraband terms appear and the result is an analytic power 
series containing only even powers of 

The interband term in Kambersky’s formula can be given a very simple interpretation in terms of Fermi’s ’’golden 
rule” for transition probability 0 . This corresponds to second order perturbation theory in the spin-orbit interaction. 
The decay of a uniform mode (q = 0) magnon into an electron-hole pair involves the transition of an electron from 
an occupied state to an unoccupied state of the same wave-vector. This is necessarily an interband transition and 
the states involved in the matrix element are unperturbed, that is calculated in the absence of SOC. A quite different 
approach has been adopted to try and find a physical interpretation of Kambersky’s intraband term ( |,i). This 
employs Kambersky’s earlier ’’breathing Fermi surface” model ( ilH) whose range of validity is uncertain. 
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We now briefly discuss the consequences of a breakdown of spatial inversion symmetry so that total orbital angular 
momentum is not quenched. In general response functions such as w) with one reversed spin are no longer 

zero to first order in Hence iSii is not given to order just by the first two terms of (l?7l) but involves further terms 
which depend explicitly on the long-range Coulomb interaction. Consequently a has a similar dependence which 
does not emerge from the torque-correlation approach. In Appendix |A] it is pointed out how the direct proof of the 
Kambersky formula breaks down in the absence of spatial inversion symmetry. 


VII. EXPERIMENTAL ASPECTS 

The inclusion of intraband terms in the Kambersky formula, despite their singular nature, has gained acceptance 
because they appear to explain a rise in intrinsic damping parameter a at low temperature which is observed in some 
systems [l5|. The calculated intraband contribution to a is proportional to the relaxation time r and it is expected 
that, due to electron-phonon scattering, r will increase as the temperature is reduced. This is in qualitative agreement 
with data [l^ for Ni and hep Co. Also a small 10% increase in a is observed in Co 2 FeAl films as the temperature is 
decreased from 300 K to 80 K . However in Fe the damping a is found to be independent of temperature down to 
4 K [I^. Very recent measurements [l^ on FePt films, with varying antisite disorder x introduced into the otherwise 
well-ordered structure, show that a increases steadily as x increases from 3 to 16%. Hence a increases monotonically 
with scattering rate 1/r as expected from the Kambersky formula in the absence of intraband terms. Furthermore 
for X = 3% it is found that a remains almost unchanged when the temperature is decreased from 200 to 20 K. Ma et 
al therefore conclude that there is no indication of an intraband term in a. From the present point of view the 
origin of the observed low temperature increase of a in Co and Ni is unclear. Further experimental work to confirm 
the results of Bhagat and Lubitz [I^ is desirable. 

The second unusual feature of the intraband term in Kambersky’s formula for a is its dependence on the SOC 
parameter This contrasts with the dependence of the interband contribution which has been observed in a 
number of alloys at room temperature [^. Recently this behaviour has been seen very precisely in FePdi_xPtx alloys 
where ^ can be varied over a wide range by varying x [isj ]. Unfortunately this work has not been extended to the 
low temperature regime where the dependence, if it exists, should be seen. It would be particularly interesting to 
see low temperature data for NiPdi_xPtx and CoPdi_xPtx since it is in Ni and Co where the intraband contribution 
has been invoked to explain the low temperature behaviour of a. From the present point of view, with the intraband 
term absent, one would expect behaviour over the whole temperature range. 


VIII. CONCLUSIONS AND OUTLOOK 


In this paper we analyse two methods which are used in the literature to calculate the damping in magnetization 
dynamics due to spin-orbit coupling. The first common approach is to employ Kambersky’s formula for the Gilbert 
damping parameter a which delivers an infinite value for a pure metal if used beyond second order in the spin-orbit 
parameter The second approach [l^ is to calculate numerically the line-width of the ferromagnetic resonance seen 
in the uniform transverse spin susceptibility. This is always found to be finite, corresponding to finite a. We resolve 
this apparent inconsistency between the two methods by an analytic treatment of the Costa-Muniz approach for the 
simplified model of a ferromagnetic metal with d-bands only. It is shown that this method leads to the Kambersky 
result correct to second order in ^ but Kambersky’s intraband scattering term, taking the non-analytic form is 

absent. Higher order terms in the present work are analytic even powers of The absence of Kambersky’s intraband 
term is the main result of this paper and it is in agreement with the conclusion that Ma et al [l^ draw from their 
experiments on FePt films. Further experimental work on the dependence of damping on electron scattering-rate and 
spin-orbit parameter in other systems is highly desirable. 

A secondary conclusion is that beyond second order in ^ some additional physics arises which has not been remarked 
on previously. This is the role of long-range Coulomb interaction which is essential for a proper treatment of the 
longitudinal susceptibility and charge response to which the transverse susceptibility is coupled by spin-orbit interac¬ 
tion. Costa and Muniz [l^ stress this coupling but fail to introduce the long-range Coulomb interaction. Generally, 
however, it seems unnecessary to go beyond second order in ^ [l3, [3 most bulk systems Kambersky’s for¬ 

mula, with electron states calculated in the absence of SOC, should be adequate. However in systems without spatial 
inversion symmetry, which include layered structures of practical importance, the Kambersky formulation may be 
inadequate even to second order in The long-range Coulomb interaction can now play a role. 

An important property of ferromagnetic systems without inversion symmetry is the Dzyaloshinskii-Moriya inter¬ 
action (DMI) which leads to an instability of the uniform ferromagnetic state with the appearance of a spiral spin 
structure or a skyrmion structure. This has been studied extensively in bulk crystals like MnSi 37| and in layered 
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structures [s^. The spiral instability appears as a singularity in the transverse susceptibility x(q, 0) at a value of q 
related to the DMI parameter. The method of this paper has been used to obtain a novel closed form expression for 
this parameter which will be reported elsewhere. 

In this paper we have analysed in some detail the transverse spin susceptibility Xj.t-i-j. but combinations of some of 
the 15 other response functions merit further study. Mixed charge-spin response arising from spin-orbit coupling is 
of particular interest for its relation to phenomena like the spin-Hall effect. 


Appendix A: A direct derivation of the Kambersky formula 


In this appendix we give a rather general derivation of the Kambersky formula for the Gilbert damping parameter 
a with an emphasis on its restriction to second order in the spin-orbit interaction parameter 
We consider a general ferromagnetic material described by the many-body Hamiltonian 

H = Hx+ H,nt + He,t (Al) 

where Hi is a one-electron Hamiltonian of the form 

Hi = Hk+ + V. (A2) 


Here Hk is the total kinetic energy, is the spin-orbit interaction, H is a potential term. Hint is the 

Coulomb interaction between electrons and Hg^t is due to an external magnetic field in the z direction. Thus 
Hext = —S^bex where bex = ^fisBex, as in (IMll . and is the z component of total spin. Both H^° and V can contain 
disorder although in this paper we consider a perfect crystal. Following the general method of Edwards and Fisher 
we use equations of motion to find that the dynamical transverse susceptibility x{^) = X_h(0;a;) satisfies [s^ 


X(w) = 


2 (^^) 

hbj - bex 




(A3) 


where 




(A4) 


with F = [S' This follows since S commutes with other terms in Hi and with Hint- For small uj, x is 

dominated by the spin wave pole at hio = bext + bSuj where ^ so that 


, ^ 2(S-) 

h(UJ — doj) — bex 

Following [ 3 ^ we compare (IA3I) and (IA5I) in the limit hSuj hw — bex to obtain 

- 2{S^)h6u; = eiXFiu;) - r\{F-,S+])) = ^^1 - lim(i([F-, S+]))] 


(A5) 


(A6) 


correct to order It is important to note that the limit ^ —S' 0 within the bracket must be taken before putting 
bu! = bex ■ If we put huj = bex fifst it is clear from (IA3I) that the quantity in brackets would vanish, giving the incorrect 
result 5uj = 0. Furthermore it may be shown [M. Cinal, private communication] that the second term in the bracket 
is real. Hence 


^{huj) = lim 3[x|, “(w)]. (A7) 

Kambersky derived this result, using the approach of Mori and Kawasaki ( [4l| El), without noting its restricted 
validity to second order in This restriction is crucial since, as discussed in the main paper, it avoids the appearance 
of singular intraband terms. Oshikawa and Affleck emphasise strongly a similar restriction in their related work on 
electron spin resonance (Appendix of El)- 

Equation (TATI) is an exact result even in the presence of disorder in the potential and spin-orbit terms of the 
Hamiltonian. In the following we assume translational symmetry. 
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To obtain the expression dm) for a = S(Sw)/&ex, which is equivalent to Kambersky’s result (l62l) . it is necessary to 
evaluate the response function in tight-binding-RPA. Using (1551) we find 


k/ii/ 


(A8) 


Hence 


xr=EE[^ /lly- G/iti^t,/3t“4-)] (-^9) 

ct/S 


where 


Gfiava' ,0TaT' (( ^ ^ *"kuo-i ^ ^ 




(AlO) 


The Green function G is to be calculated in the absence of SOC = 0). Within RPA it satisfies an equation of the 
form 


G 


^ay<7' ^proLT 


^ ^0 

' ^paua' ,l3TOi.r' 


E E 

Pl(TxVia[ P2(72^2<72 


^0 T7- 

^ paua' ,piaiUi(T[ ,/3rar' 


(All) 


where G° is the non-interacting (Hartree-Fock) Green function and 

i^l<Til'i(Tj,/i20’2l'2a'2 i^l0’J(T2 0-2 *^^*21^2 

with U(q) given by (ITTl) and (fT^ . Hence 


(A12) 


_ ^0 

Or/i(Ti/cr',/3raT' — fiaua' ,/3TaT 


- E 


E G' 

z_^ I 


rO 

aaua' 






, n 


(A13) 


The form of the interaction V given in (IA12I) is justified by the connection between (IA13I) and (fTUl) . with q = Q. To 
see this connection we note that Xaa’rr’ = J2fj,v Gfiafia' ,vtvt' and that (IA13I) then leads to (flUl) which is equivalent to 
RPA. On substituting (IA13I) into (IA9I) we see that the contributions from the second term of (IA13I) contain factors 
of the form 


E T ^ G^ \ ' r - (^0 

fiiypi 


(A14) 


We now show that such factors vanish owing to quenching of orbital angular momentum in the system without SOG 
= 0). Hence the Green functions G in (IA9I) can be replaced by the non-interacting ones G°. The non-interacting 
Green functions G° are of a similar form to x° in (l40l) and for ^ = 0 may be expressed in terms of the quantities 
fliT/i = (k/r|kn) where Ikn) is a one-electron eigenstate as introduced in section IVll Hence we find, in the same way 
that (Hil) emerged, 


E T ^ G*^ 


EE^ 

pi' kn 


/kn'1' /knj. 


/X + bea^-huj + ir]' 


(A15) 


Also by closure 


E = (kn|L^|kn) = 0, (A16) 


the last step following from quenching of total orbital angular momentum. The proof that the second expression 
in (jAldp vanishes is very similar. 

Hence we can insert the non-interacting Green functions G° in (IA9I) and straight-forward algebra, with use of (TATI) , 
leads to (l58)) . At the end of section lyTl this is shown to be equivalent to the Kambersky formula for a. We emphasize 
again that the present proof is valid only to order so that the one-electron states used to evaluate the formula 
should be calculated in the absence of SOG. 

This proof relies on the quenching of orbital angular momentum which does not occur in the absence of spatial 
inversion symmetry. When this symmetry is broken it is not difficult to see that the second term of (|A13I) gives a 
contribution to the first term on the right of (IA9l) which contains the q = 0 spin-wave pole and diverges as fiw —> bf,^. 
Hence the proof of the torque-correlation formula (IAtI) collapses. The method of section |VT] must be used as discussed 
at the end of that section. 
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Appendix B: Elements of S 


The element S'!! of the matrix S is given in (1271) . The remaining elements are given below. 

5'i2 = -Uxifif + {U/A)[{X + - (F + X|-|-4,4,)x°4,4,tX°ttt 

-{Y + + (X + x°-|-'t-|-)XuitX°tu] 


(Bl) 


-521 = -Uxnn + + Xuu)xtmXTm “ (^ + Xttu)4tnXtttt 

— (Y + x°4,tt)XtttiXt4.4.i + (-^ + Xtttt^XutfXtiu] 


(B2) 


522 — 1 — C^XtUt + (^/^)[(^ + X°4,u)XttttXt4.Tt “ (^ + Xttu)x°utXtitt 


(B3) 
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